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Strong mixing between photons and axion-like particles in the magnetic fields of clusters 
of galaxies induces a scatter in the observed luminosities of compact sources in the clus- 
ter. This is used to construct a new test for axion-like particles; applied to observations 
of active galactic nuclei it is strongly suggestive of the existence of a light axion-like particle. 



1 Introduction 

An Axion-Like Particle (ALP) is any scalar or pseudo-scalar field which couples to the kinetic 
terms of the photon. The pseudo scalar coupling to photons is identical to that of the axion; 
C D -^e^j^yXpF'-'-'^ F^P , and a scalar field couples through the Lagrangian term; C D ^F^^.i^'"'. 
The presence of contact interactions between ALPs and photons means that ALPs affect the 
propagation of photons through a magnetic field. In such an environment a photon can oscillate 
into an ALP with probability [1] 



Here z is the distance traveled, A(z) = m^^^z/^ui and tan26' = 2Buj/Mm'^j:j. rn^^^ = |?7i^ — 
cjpl, 7710 is the ALP mass, ujp the plasma frequency of the medium, u the photon frequency 
frequency, B the magnetic field strength and M the strength of the photon- ALP coupling. 

In these proceedings we describe a new test for ALPs which looks for the effects induced 
by strong ALP-photon mixing on the luminosity of astronomical objects observed through the 
magnetic fields of galaxy clusters. Our results apply to ALPs with masses < 10~^^ eV. 
The constraints on the couplings of such ALPs are: 10^^ GeV < M for pseudo-scalars [2], and 
10^^ GeV < M for scalars [3]. However a subclass of scalar ALPs known as chameleonic [1] 
ALPs avoid these constraints because their mass depends on the local density, their coupling is 
required to satisfy 10^ GeV < M [5]. 
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2 Astronomy with ALPs 



The magnetic fields of galaxy clusters fluctuate on many different scales. However, at the 
high frequencies we consider in what follows the simple cell magnetic field model can be 
shown to give the same results for ALP-photon mixing as modeling the variations in the mag- 
netic field with a power spectrum. The cell model of the magnetic field assumes the field is 
made up of a large number of equally sized magnetic domains. The magnitude of the field 
strength is the same in each domain but the orientation of the field is randomly chosen. 



When the probability of mixing between ALPs 
and photons is large the system of photons and 
ALPs can be evolved through a large number of 
randomly oriented magnetic domains analytically; 
this is known as the strong mixing limit. If L is 
the size of a magnetic domain, N the number of 
domains traversed and P = P{L) is the probabil- 
ity of photon to ALP conversion in one magnetic 
domain, we say that strong mixing occurs when 
NP > 1 and NA{L) < tt/2. In this limit the 
probability of mixing is large, and frequency in- 
dependent. Strong mixing occurs in the magnetic 
fields of galaxy clusters for x-ray or gamma-raj 
photons if M < 10^^ GeV, assuming < 
Particles with such masses and couplings are al- 
lowed by current observations for pseudo-scalar 
fields and for chameleonic scalars. 

As photon number is not conserved photon- 
ALP mixing will change the apparent luminosity 

of objects observed through the cluster. We define the attenuation factor to be the ratio of the 
flux of photons after passing through N domains to the initial flux of photons; C = Ij{N)/I^{0). 
Then in the strong mixing limit, assuming no initial flux of ALPs, the mean value for C is 
C = 2/3 [7], and its probability distribution is 8J 
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Figure 1: Simulated fingerprint for best fit 
Gaussian model 
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where a = (1 +po)/{l — po), c± — min (c, (1 ±po)/2) and po is the initial polarization of the 
photons. This probability distribution has an unusual shape, and is very asymmetric about the 
mean. In the next Section we show that this can be exploited as a new test for ALPs. 



3 Searching for ALPs with luminosity relations 

To use the shape of the probability distribution ([2]) to look for ALPs we would need to know the 
high energy photon flux for a class of astronomical sources. We do not currently know of any ob- 

^ The strength of the magnetic field is B « 1 — 10 fiG, the size of a magnetic domain is L « 1 kpc and for a 
typical source inside the cluster we expect the light observed from that source to have traversed A'^ si 100 — 1000 
magnetic domains [B]. The plasma frequency in the intracluster medium is Wp fa 10~^^ eV 
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jects that are standard candles in x- or gamma-rays, however, for certain classes of object there 
exist luminosity relations which correlate the high frequency luminosity of an object with a fea- 
ture of its low energy spectrum. At low frequencies light mixes weakly with ALPs and hence we 
assume that low energy observables are not affected by ALPs at leading order. Therefore lumi- 
nosity relations can be used to normalize the high energy flux, so that the effects of ALPs are ob- 
servable. 

The relations typically take the form 

logio Y,=a + b logio Xi + Si, (3) 

where Yi is the high energy luminosity, and Xi is 
the low energy feature of the spectrum for the i- 
th object in the survey. Si represent the scatter 
in individual measurements, it is standard in as- 
tronomy to assume they are normally distributed; 
Si = aSi where S - iV(0, 1). If the high fre- 
quency light mixes strongly with ALPs this will 
appear as an additional contribution to the scat- 
ter Si = aSi — logj^g Ci , where the Ci are described 
by the probability distribution function 

For a given data set we use the likelihood ratio 
test to see if the data prefer strong ALP-photon 
mixing, or the null hypothesis of Gaussian noise. 
We find the values of the parameters a, b and a pig^j-c 2: Simulated fingerprint for best fit 
which maximize the likelihood of each hypothe- ALP strong mixing model 
sis given the data, and then compare these two 

maximum likelihoods with the ratio r(po) = 21og(Li(po)/Lo), where Li(po) is the maximum 
likelihood allowing for strong ALP-photon mixing and Lq is the maximum likelihood for models 
where the scatter is purely Gaussian. The two hypothesis have the same number of param- 
eters and therefore r{pQ) is equivalent to the Bayesian Information Criterion. Negative r{pQ) 
is evidence against ALP strong mixing, and positive r{po) is evidence for ALP strong mixing. 
|r(po)| > 6 is considered strong evidence, |?'(po)| > 10 is considered very strong evidence. 

4 Results from active galactic nuclei 

To apply the test developed in the previous section we require a class of compact objects 
within galaxy clusters that emit x-ray or gamma-ray light and for which luminosity relations 
exist correlating the high energy luminosity with a feature of the low-energy spectrum. Active 
galactic nuclei (AGN) satisfy these requirements. For AGN a luminosity relation has been 
established between the 2 keV x-ray luminosity and the 5 eV optical luminosity. We have 
observations of 77 AGN from the COMBO-17 and ROSAT surveys and 126 objects from 
the SDSS survey [TU]. 

Applying the likelihood ratio test described in the previous section to these results we find 
r{po < 0.5) « 25, where the expectation from AGN physics is that po < 0.1 [TT] . 

As a qualitative check of this result we plot fingerprints of the data. To do this we construct 
10^ new data sets, of the same size as the original, by bootstrap re-sampling (with replace- 
ment) of the original data set. For each data set we calculate the statistical moments of the 
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distribution fcm(si) = [{l/Np) ^ ■ s™]^/™ where Si = log;^Q Yi — (a + &log;^Q Xi). These moments 
parametrize the shape of the probabihty distribution. Fingerprints of the data are then his- 
togram plots of ki vs. kj for the resampled data sets. Figures [1] and [2] show example fingerprints 
for simulated data respectively without and with the effects of strong ALP-photon mixing. 
Darker regions indicate a higher density of points. 
Figure [3] shows the same plot for the data obtained 
from observations of AGN. 

Comparing Figures [2] and [3] there is a clear 
qualitative similarity between the shapes of the 
predicted and observed distributions. This simi- 
larity persists when higher moments of the distri- 
bution arc plotted. 

The astrophysics underlying the luminosity re- 
lation for AGN is not known, and wc cannot rule 
out that a combination of standard physical pro- 
cesses in the AGN conspires to mimic the effects 
of ALP-photon mixing. It can be shown, however, 
that the scatter in the luminosity relation is not 
redshift dependent, therefore the observed scat- 
ter is not due to evolution effects or an incorrect 
choice of cosmological model. 
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Figure 3: Fingerprint from AGN data 
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